Phonon— Mediated Anomalous Dynamics of Defects 
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Dynamics of an array of line defects interacting with a background elastic medium is studied 
in the linear regime. It is shown that the inertia! coupling between the defects and the ambient 
phonons leads to an anomalous response behavior for the deformation modes of a defect-lattice, 
in the form of anisotropic and anomalous mass and elastic constants, resonant dissipation through 
excitation of phonons, and instabilities. The case of a single fluctuating line defect is also studied, 
and it is shown that it could lead to formation of shock waves in the elastic medium for sufficiently 
high frequency deformation modes. 
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Topological defects have a wide range of interest in 
many branches of physics. Vortices in superconductors, 
superfluids, and Bose-Einstein condensates, crystal de- 
fects in solid bodies, topological excitations in spin sys- 
tems with continuous symmetry, and defects in liquid 
crystals are famous examples of their ubiquity jl) . While 
they play a pivotal role in our fundamental understand- 
ing of the superfluid Q and melting || transitions in 
two dimensions, they are also important in determining 
the properties of real materials — for example they are 
responsible to a large degree for the plastic behavior of 
solid bodies ||. 

Over the past decade, the problem of defect dynam- 
ics has attracted a considerable attention mainly in de- 
scribing the dynamics of defect-mediated phase transi- 
tions [|5|-0 , and in determining the physical properties of 
vortex-lattices in type-II supercondoctors How- 
ever, the interaction of defects with the excitations of 
the elastic background, or phonons, is not taken into ac- 
count in these studies, because in the context of equi- 
librium statistical thermodynamics the defects are com- 
pletely separated from the phonons and there is no net 
coupling between the two degrees of freedom . Quan- 
tized phonons, however, are shown to interact with de- 
fects leading to interesting effects ]TT|-p^[. The interest- 
ing observation of a collection of quantized vortices in 
Bose-Einstein condensates has also raised new questions 
about the dynamical behavior of defects H] . 

Here, we consider the classical dynamics of a collec- 
tion of line defects in the form of a regular array, which 
are dynamically fluctuating in an ambient elastic medium 
away from thermodynamic equilibrium. The inertial cou- 
pling between the defects and the phonons of the back- 
ground medium is shown to lead to a plethora of novel ef- 
fects in the long time- and length-scale collective dynam- 
ics of the defect-lattice: (i) The deformation modes of 
the defect-lattice acquire an anomalous transverse mass, 
which diverges for modes that are in resonance with the 
background phonons, and a finite longitudinal mass, (ii) 
There is a corresponding anomalous transverse elastic 
constant in the direction parallel to the line defects, as 
well as a corresponding finite longitudinal clastic con- 



stant, (iii) An effective anomalous friction appears at 
the resonance, signalling the transfer of mechanical en- 
ergy from the defect-lattice to the elastic phonons. (iv) 
The defect-lattice becomes intrinsically unstable for fre- 
quencies higher than the phononic resonance frequency 
for each wavevector, due to the enhanced excitation of 
phonons. (v) The elastic moduli in the perpendicular 
direction are anomalous, and show an intrinsic instabil- 
ity for the shear modes. The case of a single fluctuating 
line defect is also studied. It is shown that deformation 
modes along the line defect with a phase velocity higher 
than the bulk velocity of phonons lead to the creation of 
shock waves. The propagation of phonons in the peri- 
odic matrix of the defects is also studied, and it is pro- 
posed that such arrangements may lead to formation of 
frequency gaps in the phononic band structure. 

We consider a simple scalar elasticity described by the 
field 9(r, t), which could be a component of the displace- 
ment field in crystals, or the phase of the order param- 
eter for superfluids. The field is then conveniently de- 
composed into two parts as 6(r, t) = # p h(r, t) + 9dei( r > t): 
(i) a singular part #def(r, t) that is a solution of the de- 
fect condition of the form <f dl ■ V0 = 27m, where the 
integral is taken over any closed path around the defect 
and n is the corresponding winding number or the topo- 
logical charge of the defect, and (ii) a nonsingular part 
# p h(r,£) that describes the phononic degrees of freedom 
in the elastic medium. While such a representation as- 
sumes a vanishing size for the defects, we should keep in 
mind that in reality they always maintain a finite core 
size, which is of the order of the atomic lattice constant 
a in crystals. With the above decomposition, we then 
consider an action as 



A= I dttfr 



<9 f 6>dcf) 2 — -^-(V6> p h 



V6> 



def J 



(1) 



where the coefficient p is a (linear) mass density, and J 
is a stiffness coefficient. The quantity c = \fjjp is the 
phase velocity of sound waves in the medium. 
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FIG. 1. Array of fluctuating line defects in an elastic 
medium. The deformation modes of the defect lattice are 
coupled to the phonons of the background medium. 

To study the many-body effects in the dynamics of 
line defects mediated by phonons, we consider a system 
consisting of many line defects (with n = 1) that are ar- 
ranged on a triangular lattice (when at equilibrium) as 
shown in Fig. 1. In this case the defect field is given as 

dcf (r ± , {Q(R, z, t)}) = are tan V Z nip ~ t ^ 

a; (aji [ti,, z, t) 

where rj_ = (x,y), and Q(R, z,i) = 
(Qi(R, z, t), Q2(R, z, t)) is the two dimensional position 
vector of the core of a line defect in the x-y plane, which is 
fluctuating around a point R of a regular two dimensional 
lattice (see Fig. 1). The corresponding Euler-Lagrange 
equation for the phonon field then reads 



[d 2 -c 2 V 2 ] 9 ph (r,t) = - [d? 



c 2 dl] def (r ± ,{Q(R, 

(3) 



and {G} represents the reciprocal lattice of the original 
defect-lattice {R}. Note that the kernel consists of two 
contributions: (i) a term which does not depend on fre- 
quency and comes from the usual long-range logarithmic 
interaction between the line defects, and (ii) a term with 
dynamical origin describing an additional interaction be- 
tween the defects mediated by phonons. 

The summation over all the reciprocal lattice vectors 
of the two dimensional lattice in the above kernel is diffi- 
cult to perform exactly. However, we can study the long 
time- and length-scale elasticity of the defect lattice by 
focusing on the small frequency and wavevector limit in 
Eq. (§). This can be achieved by separating the G = 
term and expanding the other terms in powers of q/G, 
k z /G, and lo/(cG). In this case the summation over G 
can be performed and we obtain an analytical expression 
for the kernel to the leading order as 

Xa/3(q, k s ,u) = (M T u 2 - K\k 2 z - K^q 2 ){b a!) - q a qp) 
+ (M L to 2 - K\k 2 z - Kiq 2 ) q a qp 
+ iuQt {S a /3 - q a qp), (6) 

where the effective transverse and longitudinal mass den- 
sities are given as 



Note that Vj_#dcf = 0. The above equation implies that 
the system of fluctuating defects acts as a complicated 
source for phononic waves — each defect can transmit or 
absorb phonons from the medium during its motion. We 
will see below that the net result of this coupling is a 
dissipative effect in which the moving defects lose their 
energy in the form of phonon radiation. 

To study the collective dynamics of the defects in the 
presence of phonons, we should eliminate the phononic 
degrees of freedom in the action, and obtain an effective 
action for the defect degrees of freedom {Q(R, z,t)}. For 
this purpose, we assume that the line defects are only 
slightly distorted from their equilibrium configurations, 
so that their positions can be described as Q(R, z,t) = 
R+u(R, z, t) with u(R, z, t) being much less than the de- 
fect lattice constant b. We then solve the wave equation 
for the phonons [Eq. (||)] perturbatively in the deforma- 
tion field u(R, z, t), and obtain the effective action up to 
the leading order as 



M T (q, k z ,uj) 



q 2 + k 2 — (uj/c) 2 K 



P M L =? (7) 



the effective transverse and longitudinal elastic moduli in 
the parallel direction as 



K { ^(q, k z ,u>) 



J 



q 2 + k 2 — (lo/c) 2 k 



7 T ^> K l - ^2> ( 8 ) 



the effective transverse and longitudinal elastic moduli in 
the perpendicular direction as 



(9) 



and, finally, the effective transverse friction coefficient as 



C,r{q,k z ,uj) 



ttJc q 



2„1 



S(u 2 - c 2 k 2 - c 2 q 2 ). (10) 
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In the above, parallel and perpendicular are denned with 
respect to the line defects, and the screening length k^ 1 
is defined as 



-y- 



n\n(b/a) 
G* 2 



(11) 



where G* is the lattice constant of the reciprocal lattice, 
which for a triangular lattice of line defects is given as 
G* = ^=£. Note that the real (imaginary) part of the 

above kernel is even (odd) in u> and this ensures that the 
response is causal. 

The above results show that the dynamics of a defect 
lattice in an elastic medium is anomalous. While the 
longitudinal deformation modes acquire finite mass and 
elastic modulus in the parallel direction due to coupling 
with phonons, the corresponding mass and elastic modu- 
lus in the parallel direction for the transverse modes are 
frequency and wavevector dependent, and diverge in the 
small wavevector and frequency limit. Moreover, whereas 
the transverse mass and elastic modulus in the parallel 
direction are positive for lo 2 < c 2 (q 2 + k 2 ), they diverge 
when the frequency and wavevectors of the deformation 
modes hit that of a possible phononic excitation in the 
background. In this case the defect-lattice dynamics be- 
comes dissipative — as manifested by the appearance of 
an imaginary part in the response kernel that defines an 
effective friction coefficient — which corresponds to trans- 
fer of mechanical energy from the defect lattice to the 
phonons. For lo 2 > c 2 (q 2 + k 2 z ) 1 the transverse mass and 
elastic modulus in the parallel direction become nega- 
tive, signalling an instability in the defect-lattice due to 
the resonant coupling with the background phonons. The 
transverse and longitudinal elastic moduli in the perpen- 
dicular direction are also wavevector dependent, and they 
reveal an inherent instability in the deformation modes 
that is related to the well-known instability of a system 
of charges in electrostatics; often alluded to as the Earn- 
shaw's theorem P,^5[. 

It is also interesting to consider the case where all the 
line defects are frozen (to straight lines) and only a single 
defect is fluctuating, which can be achieved by assuming 
that u a (q, k z , lo) does not depend on q. We can then 
combine the restricted q integration in the first Brillouin 
zone and the summation over the reciprocal lattice vec- 
tors G to a free integration over q, and calculate the 
response kernel Xsi{k z ,Lo)S a/3 = J B Z Xa/3(q, k z ,oS). 

This yields 



Xs\(k z ,uj) = p s \uj 2 - K sl k 2 + iujQ 



(12) 



with the effective linear mass density for a single line 
defect given as 



p s i(k z ,uj) = - p In 



2 r "~ V \uj 2 - c 2 k 2 
the effective elastic modulus as 



(13) 



K sl (k z ,uj) 



J In 



c 2 k 2 \ 



and the effective friction coefficient as 
- 2 (lo 2 - c 2 k 2 ) 



(sl{k z ,w) 



7T 



P 



ho 



e(^-c^) 



(14) 



(15) 



In the above, the Debye frequency lo-q — c/a is a high fre- 
quency cutoff in the system, and represents the Heav- 
iside step function. 

The effective mass density and elastic modulus for the 
line defect are logarithmically larger than the nominal 
mass density and stiffness of the material in the core of 
the vortex, and they are also frequency and wavevector 
dependent. While they are positive for all frequency and 
wavevectors, they logarithmically diverge when lo = ck z , 
and in particular in the limit lo — > and k z — > 0. 

In the case where the phase velocity of waves on the 
line defect Lo/k z is greater than the phonon velocity in 
the system c, the response kernel in Eq. (12) develops an 
imaginary part, which is a reflection of a transfer of me- 
chanical energy from the fluctuating defect to the elastic 
medium in the form of radiating phonons. This effect 
is similar to the Cherenkov radiation, where a charged 
particle moving in a dielectric medium with superlumi- 
nal velocity (i.e. a velocity greater than the velocity of 
photons in that medium) radiates electromagnetic waves. 

The time-averaged energy dissipation rate for the fluc- 
tuating line defect can be calculated as 



P 



dk z dio 
~2tt2it 



( s i(k z ,to) lo 2 \u(k z ,Lo)\ 2 



(16) 



To explicitly check that the dissipation is due to phonon 
radiation, we calculate the rate of change in energy for 
the phonon field as 



dE, 



dt 



d 3 r -|t 00 , 
at 



(17) 



using the 00-component of the stress tensor |T(| that 
can be calculated from the phonon field action in Eq. (Q). 
The result comes out exactly equal to the dissipation rate 
as given in Eq. ( |l6| ) above. 

To estimate the radiated power, we consider an ex- 
ample in which the defect is undergoing a solid harmonic 
motion described via u(k z , lo) — u (27i ) 2 S(k z ) [5 (lo+lo ) + 
5(lo — loq)]/2. The time-averaged radiated power per unit 
length from this oscillating line defect is then calculated 
from Eq. © as P/L = 7r 2 pugw^/4. 

Finally, we note that the periodic arrangement of the 
defect lines on a lattice changes the spectrum of phonons, 
as in the case of electrons in metals. This is manifested in 
the singular points (poles) in the dynamical term in the 
response kernel of Eq. dq) whose structure determines 
the dispersion relation for the phononic excitations. The 
dispersion relation f2(q) = \Juo(q) 2 — c 2 k 2 — c|q — G| 
should respect the symmetry of the reciprocal lattice G. 
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As an example, the phonemic dispersion diagram for the 
case of a triangular lattice of line defects is plotted in 
Fig. 2. The unit vectors for the triangular lattice are 
6(1,0) and and the unit vectors for the corre- 

sponding reciprocal lattice are G*(^, — A) and G*(0, 1), 
where G* — 4=r is the lattice constant of the reciprocal 




FIG. 2. Di agram of pho nemic dispersion relation defined 
via fi(q) = ■J 'u>(q) 2 — c 2 k 2 = c|q — G in a triangular lattice 
of line defects, with flo = cG* — ^=|. The inset shows the 
first Brillouin zone of the triangular lattice. 

While the band diagram in Fig. 2 has no gaps, one can 
argue that an interaction of the phonons with the inte- 
rior of the defects that changes the propagation velocity 
of phonons from c to c/y/ e(r) in the vicinity of the de- 
fects H can produce a gap in the band diagram, similar 
to the case of photonic crystals (Tt]]. A simple calcula- 
tion (similar to the estimation of the gap for electrons 
in a weak periodic potential B ) then shows that such an 
interaction can lead to frequency gaps in the zone bound- 
ary of the order of A ~ cG* (a/b), where a is the core size 
of the defect |Q. This may suggest that a rather dense 
array of defects may act as a "phononic crystal," which 
could be of potential interest for applications. 

In conclusion, we have shown that the dynamical in- 
teraction of defects with phonons can lead to renormal- 
ization of mass and dissipation, which is analogous to 
the radiation-reaction force for electrons and their mass 
renormalization due to their interaction with electromag- 
netic fields @. 
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